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ABSTRACT 

The hot Jupiter WASP-12b is a heavily irradiated exoplanet in a short pe- 
riod orbit around a GO-star with twice the metaUicity of the Sun. A recent 
thermochemical equihbrium analysis based on Spitzer and ground-based infrared 
observations suggests that the presence of CH4 in its atmosphere and the lack of 
H2O features can only be explained if the carbon-to-oxygen ratio in the planet's 
atmosphere is much greater than the solar ratio ([C]/[0] = 0.54). Here, we use a 
1-D photochemical model to study the effect of disequilibrium chemistry on the 
observed abundances of H2O, CO, CO2 and CH4 in the WASP-12b atmosphere. 
We consider two cases: one with solar [C]/[0] and another with [C]/[0] = 1.08. 
The solar case predicts that H2O and CO are more abundant than CO2 and 
CH4, as expected, whereas the high [C]/[0] model shows that CO, C2H2 and 
HCN are more abundant. This indicates that the extra carbon from the high 
[C]/[0] model is in hydrocarbon species. H2O photolysis is the dominant dise- 
quilibrium mechanism that alters the chemistry at higher altitudes in the solar 
[C]/[0] case, whereas photodissociation of C2H2 and HCN is significant in the 
super-solar case. Furthermore, our analysis indicates that 2h2 is the major ab- 
sorber in the atmosphere of WASP-12b and the absorption features detected near 
1.6 and 8 micron may be arising from C2H2 rather than CH4. The Hubble Space 
Telescope's WFC3 can resolve this discrepancy, as '2h2 has absorption between 
1.51 — 1.54 microns, while CH4 does not. 

Subject headings: stars: planetary systems 



1. Introduction 



The discovery of the first transiting planet, HD 209458b (jCharbonneau et al.l 12000 



Henry et al.l I2OOOI ). opened up a new window to observe and study extrasolar planetary 
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systems. By combining transit data with rad ial velocity measu r ements, one can dete r- 
mine the mass and radius of a transit i ng pl anet flMazeh et al.ll2000l : iLaughlin et al.ll2005al Jb: 
Holman fc Murray I l2005l : lAgol et al.l |2005[ ). Apart from these physical properties of the 
planet, it has also been shown that the transmission and emission spectra from ground- and 
space-based observations can be used to place cons t raints on the atmospheric composition 



( Charbonneau et al.l 



2008 



2002 



Vidal-Madiar et al.ll2003 



Swain et al.ll2009al Jbl). brightness temperatu re (iDeming et al.ll20p5l : ICharbonneau et al. 



2004 



Tinetti et al. 



2007 



Snellen et al. 



2005 ) and even day-night temperature contrast ( Knutson et al. 2(3071) of transiting planets. 
To date, most of the planets discovered are "hot Jupit ers" flCoUier-Cameron 1120021), but re- 
cently terrestrial mass planets are also being discovered ( Leger et ahlboOQ : Charbonneau et al 
2009h . 



The observational determination of chemical species that exist in the outer atmospheric 
layers of transiting planets provides us with an opportunity to investigate the underlying 
chemistry. T ypically, thermochernical equilibrium has been assumed in rnodels of exoplanet 
atmospheres ( Burrows et al.lll997l : iFortney et al.ll2005t ISeager et al.ll2005l : iMarley et al.ll2007l : 
Madhusudhan et al .1120 111 ). This is a good assumption at the high temperatures and high 
pressures prevailing in the lower layers of close-in gas giants. Disequilibrium caused either by 
the UV flux of the host star (photochemistry) or b y eddy and molecular diffusion (vertical 



transport) have been con si dered in some models (ILiang et a. 



2006: Zahnle et al. 2009a b; Line et al. 2010: Moses et al 



20031 : ICooper fc Showman 



201ll ). These studies showed that 



disequilibrium mechanisms can be significant in determining the chemical composition of hot 
Jupiters. 

Here, we investigate the significance of photochemistry and transport in determining 
the abundances of major species, HoO, CO and CH4, observe d in the dayside thermal emis- 
sion spectru m (IMadhusudhan et al.l 1201 ll : ICroll et al.l 1201 ll ) of the transiting hot- Jupiter 
WASP-12b JHebb et al.l 120091 ). At the time of its discovery WASP-12b was the most highly 
irradiated exoplanet ( T> 25 00 K) with the largest radiu s (1.79-Rj) and the shortest orbital 
period (1.09 days). Recently, IMadhusudhan et al.l ( 1201 ll ) reported that Spitzer Space Tele- 
scope observations show strong absorption features of CH4 in the 3.6 /xm channel and CO in 
the 4.5 /xm channel, whereas weaker features were observed in the 5.8 /im channel where H2O 
absorbs. This suggests that CH4 and CO are dominant and that H2O is less abundant in the 
atmosphere of WASP- 12b. Assuming equilibrium chemistry and solar [C]/[0] = 0.54, H2O 
and CO should be the dominant species a nd CH4 and CO2 should be the least abundant. 
Therefore, a solar [C]/[0] ratio is ruled out. IMadhusudhan et al.l (120111 ) conclude that to ex- 
plain the observed abundance of CH4 and CO, WASP-12b must have [C]/[0] > 1, implying 
that it is a carbon-rich planet. 
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WASP- 12b is one of the most highly irradiated known exoplanets, so photochemistry 
could play an important role in determining its atmospheric composition. Assuming [C]/[0] = 
1, analysis of observations using equilibrium chemistry models suggests mixing ratios (with 
respect to molecular hydrogen) less than 10"''^ for H2O, greater than 10~^ for CO, ~ 10"^ 
for CH4 and less than 10~^ for CO2. Our goal in this study is to examine how the vertical 
distribution and abundances of these species are affected by photochemistry. Specifically, 
we wish to determine whether photochemical models make qualitatively different predictions 
from those of thermochemical equilibrium models. We consider two cases, one with solar 
[C]/[0] =0.54 and another with [C]/[0] = 1.08, both of which have also been studied with 
equihbrium models. 



2. Model description 

We use a one-dime nsional photoc h emical model initia l ly developed t o study primitive 



f Kastinp 


1982. 


1983; 


Zahnle 


1986: 


Kastinp 


1990) 



been modified to suit the hot Jupiter temperature regime by including "backwards" chemical 



Line et al. 


2010; 


Moses et al. 


2011) 



pled ordinary differential equations (ODEs) for the mixing ratios of all species at all heights 
using the reverse Euler method. The method is first order in time and uses second-order 
centered finite differences in space. We include the following 31 chemical species involved 
in 230 reactions: O, O2, H2O, H, OH, CO2, CO, HCO, H2CO, CH4, CH3, CHsCCHgOH, 
CH, CH2, H2COH, C, C2, C2H, C2H2, N, N2, NH, NH2, NH3, CN, HCN, H2, He, 0(^0) and 
^CH2. These species are divided into long-lived species (from O to H2), short-lived species ( 
0(^D) and ^CH2) and "inert" species (He). Both chemistry and vertical transport by eddy 
diffusion are considered for long-lived species, whereas transport is neglected for short-lived 
species. Constant mixing ratios with altitud e are as s umed for " inert" species. The reaction 
list and rate constants were obtained from IZahnld (1201 ll ) and are listed in Table [H W e 



have t aken only the "forward" reactions and corresponding rate coefficients, kf, from 



Zahnle 



(120 111 ): The reverse rate coefficients, kr, at each temperature level (grid) were calculated 
assuming thermodynamic-equilibrium: kr = kf/keq, where k^q is the equilibrium constant 
for the reaction and is given by k^q = e^'^"^^'^. Here AG° is the change in the Gibbs free 
energy for the reaction and is calculated from the Gibbs free energy of formation of reactants 
and products obtained from NIST-JANAF thermochemical tables when availabl^ (and from 



f 



http : //webbook . nist . gov/ chemistryT] 
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NASA thermobuild websit^when not available): AG° = AG°j{reactants)—AG°j (products). 
However, it should be noted that one cannot simply calculate kr as discussed above. This 
is because the rate coefficients are given in units applicable to number densities whereas the 
thermodynamic quantities (enthalpy, entropy & Gibbs energies) are calculated at a reference 
pressure (usually 1 bar). For reactions that have different number of reactants and products, 
proper pressure terms must be added to obtain accurate rate coefficients. We have appropri- 
ately include d these terms in ou r rever se rates. Discussion of these correction terms is given 
section 2.2 of Ivisscher fc MosesI J201lh . 



As lower boundary conditions, we fix the mixing ratios of the species at thermodynamic 
equilibrium values. Constant (zero) deposition velocities are assumed for the other species. 
The upper boundary condition is set to zero flux for all the long-lived spec ies. A mor e 



detailed description of the numerical scheme employed in this model is given in lPavlov et al. 
J200lh . 



The vertical grid has 100 altitude levels, ranging from km (lower boundary) to 12, 800 
km (upper boundary) in 128-km increments. The lower boundary pressure is set at 1 bar, 
and the upper boundary is fixed at 10^^ bar. Going to higher pressures is unnecessary, be- 
cause the species proffies are already close to thermodynamic equilibriu m well above the 1-bar 



l evel. For the temperature profile, we use one of the best-fit models from lMadhusudhan et al. 
(120 111 ) which has no inversion (purple curve in their Fig. 1). The pressure proffie was re- 



calculated from this temperature profile by assuming hydrostatic equilibrium and using the 
calculated mean molecular weights from the photochemical model. Vertical transport is pa- 
rameterized as eddy diffusion, as is common in one-dimensiona l photochemical models. The 
dayside eddy diffusion proffie from Fig. 1 of iLine et al.l (120101). which is o riginally obtained 
from the vertical winds from HD 189733b GCM of IShowman et al.l (120091 ). is adopted. We 
have also performed sensitivity tests by varying eddy proffies, as discussed in §11 Both the 
temperature and eddy profiles in our photochemical model are shown in Fig. [H 

For comparative purposes, we also calculate the thermodynamic equilibrium mixing ra- 
tios for all the species in the photochemical model at each altitude by solving simultaneously 
a system of chemical equilibrium equations. These equations require the total elemental 
abundances of carbon, oxygen, hydrogen and nitrogen (as we consider only compounds from 
these eleme nts) and Gibbs free eri ergies as a function of temperature. Solar elemental abun- 
dances from lAsplund et al.l (120051 ) are assumed to be our base values, but we report results 
for both solar [C]/[0] and 2x solar [C]/[0]. To calculate Gibbs free energies, the enthalpy 
of formation at the reference temperature (298 K) and entropy are needed. We then use the 



' http : //www . grc . nasa . gov/WWW/ CEAWeb/ceaThermoBuild . htm 
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expressions given in IChasd (119981 ) (p. 16) to calculate Gibbs free energy of formation for 
each species. 



We initially tested our model by attempting to reproduce the results of iLine et al. 



(120 lOl ). for the hot Jupiter planet HD189733b. The dayside temperature and eddy diffusion 
profiles were taken from their Fig. 1. The lower boundary pressure was fixed at 10 bar. 
Both the thermochemical equilibrium and photochemical model results are shown i n Fig . 
151 of appendix [K\ and a re in good agreement with the similar analyses of iLine et al.l (120101 ) 
and iMoses et al.l ( 20111 ). Moreover, our model maintains equilibrium concentrations for all 
the major species in the deeper levels (~ 10 bar), as it should at high temperatures and 
pressures. 



The star WASP-12 is a GO sta 



i§| with an effective temperature of 6350 K and twice the 
solar metallicity (IHebb et al.ll2009l ). To si mulate its spe ctrum, we used a GOV star spectrum 
from Pickles stellar spectral flux library (Pickles 1998 1^. normalized to a solar flux of 1360 
W m~^ (the value at Earth's orbit today). We then multiplied the flux at each wavelength by 
a value consistent with inverse square law of the distance to get the correct flux for WASP- 
12b. The Pickles spectra are normalized to 1 at 5556A. The fluxes from this dimensionless 
mod el spectrum are converted to W m~^nm~^ by multiplying the following flux expression 
fromlCraJ Jl992h : 



logFf 



5556 



-OAOV - 8.449 



where -F^55fi is the flux at 555 6 A and V is the visual magnitude of the star. For WASP-12, 
V = 11.69 JHebb et al.l 120091 ). In Fig. El we show the GOV star spectrum along with F2V 
and the Sun. A fixed stellar zenith angle of 50° is assum ed in all our models , the same value 
that is used in our models of Ea rth's atmosphere ( e.g. jPavlov et al.l (120011 )). This value is 



close to the value of 48° used by iMoses et al.l (120111 ) to reproduce secondary transit spectra 
in the atmosphere of HD 189733b. 



3. Results 

We consider two different cases. In the first, we assume [C]/[0] = 0.54 (solar), and 
in the second we assume [C]/[0] = 1.08 (twice solar). Fig. [3] shows mixing ratio profiles 
of some of the major species in our model, plotted against pressure for solar [C]/[0] abun- 
dance. The lower boundary in both the models is kept at 1 bar pressure (T = 2841 K) as 



• ^http : / /www ■ superwasp . org7wasp_plgLnets . htm 
''http://cdsarc.u-strasbg.fr/viz-bin/Cat7J/PASP/110/863 
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the o bserved spectral features mostly arise from pressures equal to or less than the 1-ba r 
level JPortnev et aliboosi Ixinetti et aPboOTt ISwain et al]l2009al : iMadhusudhan et alj[201ll ) . 
Dashed lines represent the profile obtained from equilibrium chemistry, solid lines from our 
photochemical model and filled squares represent mixing-ratios of respective species at the 
lower boundary. 

In the case of solar [C]/[0] (Fig. [3]), most of the oxygen and carbon is in H2O and CO. 
The chemical loss time scale (rchem) for H2O, CO and CO2 is smaller than the transport time 
scale (rtrans) ; hence, as altitude incre ases the abundances stay at their equilibrium values until 
~ 10~^ bar (jPrinn fc Barshaylll977l ). Below this pressure level (i.e., at higher altitudes) H2O 
gets photolyzed. H2O photolyzes at lower altitudes than do es CO because th e disso ciation 
energy for H2O (5.17 eV) is lower than that of CO (11.14 eV) dVung fc DeMorel Jl999h . Table 
2.4). To break this strong C-0 bond, photons of wavelength < 111.3 nm are needed (H2O 
needs photons of wavelength <239.8 to break its bond). CO photolysis is not simulated in 
our photochem ical model and it can be a source of carbon and oxygen photochemistry at 
high altitudes (ILine et al.ll2010l : iMoses et al.ll201ll ). but it should be relatively slow because 
of the small number of photons at these short wavelengths. By contrast, the photon flux is 
quite high at the longer wavelengths that can photolyze H2O (green curve in Fig. [2] inset). 

In the case of CO2, the equilibrium abundance is maintained until 10~^ bar and is set 
by the following kinetic reactions that transfer oxygen from CO and H2O to CO2: 



H2O + H ^ H2 + OH 
CO + OH ^ CO2 + H 



(2) 
(3) 



At altitudes above 10~^ bar, H2O photolysis becomes the dominant source of OH production. 
The OH then combines with CO through Eq. to produce excess CO2 (local maximum of 
solid light-blue curve > 10~^ bar). Above this level, CO2 becomes less abundant because it 
is itself photolyzed. 

The shape of the equilibrium profile for CH4 can be understood from the following 
reaction 



CO + 3H2 O CH4 + H2O 
and the corresponding equilibrium constant: 



eq 



PCH4 PH2O 
PCO Pm 



(4) 



(5) 



where 'pj' represents the partial pressure of species i. The partial pressure is related to 
the total pressure and volume mixing ratio as Pi = fi ■ P- As pressure increases (going 
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downward) from 10~^ bar, the denominator term on the right-hand side of Eq.([5]) increases. 
Temperature is constant in this region (see Fig. [1]), as is the H2O/CO ratio. Thus, in order 
to maintain equihbrium CH4 must increase with depth. Below 10~^ bar, the temperature 
starts to increase with depth. CH4 is more stable at lower temperatures and is also more 
sensitive to temperature changes than other species. Hence, it becomes less abundant in 
the 10~^ — 10~^ bar regime. At pressures above 10~^ bar, the temperature again remains 
constant, so CH4 must again increase with depth as it does in the upper atmosphere. 

The CH4 profile from the photochemical model (solid magenta curve in Fig. E]) follows 
the equilibrium profile at pressures up to ~ 10~^ bar. Above this level, Tchem ~ Ttrans (quench 
level) and CH4 remains well mixed near its equilibrium value of 10~^°. Photolysis of CH4 
occurs above 10"'^ bar (see FigJT]) mainly through the following reactions: 

CH4 + hz/ ^ CH2 + 2H, (6) 
CH4 + hz/ ^ CH3 + H, (7) 
CH4 + hz/ ^ ^CH2 + H2 (8) 

Although H is produced through CH4 photolysis, it is not enough to explain the increase 
in H abundance between 10"^ — 10"^ bar (light-blue solid curve in Figj3]). This increase in 
H is mainly due to the production of OH through H2O photolysis at this level, which then 
combines with the most abundant molecule in this atmosphere, H2, through the reverse of 
reaction ([2]). 

H20 + hz/ ^ H + OH, (9) 
H2 + OH -> H2O + H (10) 



As the OH abundance increases, more H2 is consumed and it's mixing ratio decreases 
above 10~^ bar (solid black curve in Fig. [3]). Eventually, H2O itself becomes depleted by 
photolysis, so the production of OH radical diminishes. At this point, H2 asymptotically 
reaches a mixing ratio of 10~^. The increase in H also affects the atomic oxygen abundance 
(black solid curve in Fig. [3]) by the following reaction: 

CO + H^O + CH (11) 

Note that CO photolysis, which is not included in our model, may dominate the above 
reaction in producing atomic oxygen. 

Our analysis shows that, in the [C]/[0] = 0.54 case, the abundances of major species 
(Fig. [3]) in WASP-12b's atmosphere are mainly determined by thermochemical equilibrium, 
with departures at high altitudes due to disequilibrium chemistry driven by H2O photolysis. 



- 8 - 



This is not surprising considering that H2O is far more abundant than CH4. The photolysis 
rates of H2O and CH4 as a function of pressure (altitude) for solar (blue) and super-solar 
(red) [C]/[0] are shown in Fig. [71 In the solar case, at any given height, H2O is more rapidly 
photodissociated than is CH4, as it is more abundant. In the super-solar [C]/[0] case (Fig. 
H]), 2h2 is more abundant than either H2O or CH4, and so it is photolyzed more rapidly at 
high altitudes. 



Madhusudhan et al.l (120111 ) report that the spectrum obtained from the dayside multi- 
wavelength photometry of WASP-12b is best explained if one assumes [C]/[0] > 1, using 
chemical equilibrium models. Under this assumption, the atmosphere is depleted in H2O, 
enhanced in CH4, and rich in CO. These equilibrium model profiles (dashed lines), along 
with our photochemical model results, for [C]/[0] = 1.08 are shown in Fig. HI In contrast to 
solar [C]/[0] model, most of the oxygen is now in CO (blue solid line), and CH4 (magenta 
curve) is more abundant than H2O (red curve). The switchover from an atmosphere where 
H2O and CO are the dominant species to one in which CH4 and CO become abundant 
happens precisely at [C]/[0] = 1. This transition is illustrated in Fig. [61 

In Fig. mthe abundance of H2O follows the equilibrium profile at pressures up to about 
10~^ bar in the [C]/[0] = 1.08 model. At that point, transport by eddy diffusion becomes 
faster than the chemical reaction timescale, so the equilibrium value is maintained until 
~ 10~^ bar. Photolysis begins above this level, and the abundance of H2O decreases. The 
behavior of CH4 is similar to the solar [C]/[0] case, though it is relatively more abundant in 
this high [C]/[0] model. Note that for CH4 and H2O, the photochemical mixing ratios are 
not exactly equal to the equilibrium values below 0.1 bar. The reason is as follows: At high 
altitudes (above 10~^ bar), atomic hydrogen (green solid curve in Fig. becomes a dominant 
species (more than H2). Our photochemical model uses a minor constituent approx imation 



for the diffusion coefficient in a binary mixture (Eq. (15.29). [Banks fc KockartsI (119731 )). which 
clearly is not applicable to H at this level. Due to this approximation, the mixing ratio of 
H exceeds unity above 10~^ bar, which is unphysical. Therefore, we renormalize the mixing 
ratios of each species in our photochemical model so that they sum to unity, and hence 
the equilibrium and photochemical profiles deviate slightly in the lower atmosphere. Note 
that this should not affect our conclusions in any way regarding which species are dominant 
(discussed in the next paragraph) in Wasp-12b's high [C]/[0] model. 



As can be seen from Fig. HI and also pointed by [Moses et al.[ (l2011bl ). the dominant 
hydrogen species (apart from H and H2) in this model are HCN and 2h2. Therefore, the 
photolysis of these two species drive the disequilibrium chemistry in the upper atmosphere. 
For example, in the solar model, the catalytic H2 destruction mechanism initiated by H2O 
photolysis (Eqs. dHJ & (ITOll ) was used to explain the increase in H abundance shown in Fig. 
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|3] (green solid curve) . A similar increase of H at high altitudes can be noticed in the high 
[C]/[0] case. Reactions ([9]) and ffTOj) require OH production through H2O photolysis, which 
is negligible in the high [C]/[0] model. Instead, the following reactions are important: 



2h2 + hu ^ C2H + H, (12) 

C2H + H2 ^ 2h2 + H, (13) 

HCN + hz/ ^ CN + H, (14) 

CN + H2 ^ HCN + H (15) 

The results of these reactions can be seen in Fig. HI At altitudes above ~ 10^^ bar, the 
photolysis of 2h2 and HCN produce C2H, CN and H through the above reactions. An increase 
in H can be seen as a result. The abundances of CN & C2H are not large enough below 10~^ 
bar to have a significant effect on the mixing ratio of H2. Above this level, reactions (1131) & 
f|T5|) result in the decrease of H2 mixing ratio (solid black line in Fig. H]) and corresponding 
increase of H. Further up, 2h2 and HCN become scarce and the production of C2H and CN 
diminishes, which in turn reduces the rate of production of H. Therefore, H assumes a nearly 
constant mixing ratio thereafter. 



Based on their thermodynamic equilibrium calculations, iMadhusudhan et al.l (120111 ) 
concluded that Wasp- 12b is abundant in CH4 and deficient in H2O. Our analysis indi- 
cates that both the equilibrium and photochemical models predict 2h2 and HCN are more 
abundant than CH4. Also, 2h2 has strong absorption in the range 2.98 — 3.1 microns and 
also between 7.2 — 7.9 microns, whereas CH4 has absorption features between 3.2 — 3.45 
microns and 7.3 — 8 microns. The short wavelength range for 2h2 has little overlap with the 
Spitzer 3.6 micron channerl but the longer wavelength range for both species overlaps with 
Spitzer's 8 micron channeo. In order to determine which is the dominant absorber, we have 
calculated the optical depths of the 2h2 7.5 micron band and the CH4 7.7 micron band as a 
function of pressure, as shown in FigjSl Approximate band-averaged absorption coefficients 
for these features are 2 x 10~^^ cm^ and 4 x 10~^^ cm^, respectiveljjl|. The column depths 
are taken from our photochemical model. Clearly, 2^2 has a larger optical depth than CH4 
and is the dominant absorber. Note that, Madhusudhan et "all (j201ll ) point out that 0.01 — 1 



bar pressure levels contributes most to the observed spectrum and that 2h2 is considerably 



■'Band pass range from 3.08—4.01 micron: http : / /irsa. ipac . caltech. edu/data/SPITZER/docs/irac/calibrationf il 
^Band pass range from 6.15—10.49 microns: |http: / / irsa. ipac . caltech. edu/data/SPITZER/ docs/ irac/ calibrationf : 



http : / / vpl . astro . Washington . edu/ spectra/ c2h2pnnlimagesinicrons . htm 
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more abundant than CH4 within this pressure range (Fig. H]). Therefore, future analysis of 
observations of carbon rich planets (including further analysis of WASP-12b) should consider 
higher-hydrocarbon species. 



4. Discussion 

Our analysis confirms the previous thermodjTiamic equilibrium result that a [C]/[0] > 1 
is needed to explain the observed o verabundance of CH4 in the atmosphere of WASP-12b. 



A similar conclusion was reached by iLine et al.l (120 lOf ) but for a different hot Jupiter planet 



HD 189733b. These authors varied [C]/[0] from 0.1 to 10 times the solar value, while keeping 
the total metallicity at the solar value (Fig. 6 in their paper), and examined the effect on 
thermochemcal equilibrium mixing ratios at the lower boundary. As [C]/[0] increases, most 
of the carbon in their model is in CO and CH4. At [C]/[0] = 1, H2O and CH4 switch their 
profiles just as they do in our equilibrium models of WASP-12b (Fig. E]). Although our 
equilibrium models agree qualitatively with theirs, the respective mixing ratios of the major 
species differ because of different elemental abundances and overall hotter temperatures 
(their ~ 1500 K versus our 2800 K) 

We have also performed a sensitivity test to eddy diffusion varying by 3 orders of mag- 
nitude above and below our eddy profile. For the larger case the species concentrations are 
well mixed over much of the atmosphere, deviating from the equilibrium even at relatively 
low altitudes (~ 0.1 bar). Consequently, the photolysis of 2h2 and HCN is not effective in 
producing atomic hydrogen (as mixing dominates photolysis even at h igh altitudes). On 



the ot her hand, if the eddy diffusion coefficient is small (as proposed by lYoudin &: Mitchell 



( I2OIOI )). mixing is not effective and photochemistry becomes important at mid altitudes 
(10~^ — 10~^ bar). Therefore, significant deviations from equilibrium occur at all altitudes 
above this level. 



4.1. A possible mechanism for the origin of excess carbon in WASP-12b 



The high [C]/[0] rati o in WASP-12b is une xpected, considering that the host star has 
a solar [C]/[0] ratio fsee iFossati et al.l ( 1201 Obi ) Table 2). In the standard core accretion 
model ( jPoUack et al.lll996l ). volatiles s uch as carbon an d oxygen are ex pected to remain 
unfractionated in forming giant planets (jOwen et al.lll999[ ). iLoddersI ( 120041 ) pointed out that 
Galileo probe measurements of Jupiter's atmosphere show an enriched carbon abundance of 
1.7 times solar and a depletion of oxygen by a factor of 4 (but see further discussion below). 
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To explain this result, iLoddersI fl2004[ ) proposed a model in which carbonaceous matter began 
to condense in the solar nebula beyond 5 AU, thereby providing the increased mass density 
needed for rapid core growth. By contrast, in the standard accretio n model, Jupite r forms 
just beyond the "ice line" where water ice begins to condense. In the iLodders model, 
the ice line would have been farther out, beyond the orbit of Jupiter, and this would explain 
Jupiter's apparent deficiency in O relative to C. A similar mechanism might then account 
for the high [C]/[0] ratio in WASP 12b. 



Although the iLoddersI (120041 ) model could be correct, we suspect that Jupiter formed 
beyond the ice line, with a solar [C]/[0] ratio, and that other factors are responsible for 
observed [C]/[0] enrichments in exopla nets. The Galileo Probe is widely thought to have 
descended into an infrared "hot spot" flAtreya et al.lll999l ). that is, an area of downwelling 
air that had been depleted in H2O during its uplift from below. In support of this idea, 
the H2O mixing ratio was observed to gradually i ncrease with depth d own to 20 bars, at 
which point the probe lost contact with Earth (see lAtreya et al.l (119991 ) Table 1). Further- 
more, tlujii^^erstoiins^^ deeper than 4 — 5 
bars (IGierasch et al.ll2000l : llngersoll et al.ll2000l : lAtreya et al.ll2005l ). which is consistent with 
equilibrium cloud condensatation models which predict that water clouds can form in this 
pressure range if the oxygen abundance is at least solar ( Atreya et al. 200^ 



As an alternative to the ILoddersI (120041 ) model, we suggest that the high [C]/[0] for 
WASP-12b arose because the primordial disk was depleted in oxygen abundance during the 
giant planet's migration. The carbon compounds (CH4, CO) may have been trapped in ices 
in the form of planetesimals and then accreted onto the envelope of the gas giant, resulting 
in the observed enhancement of [C]/[0] > 1. Assuming that the disk started with solar 
elemental abundances of carbon (2.26 x 10~^) and oxygen (4.20 x lO^'*), in order to obtain 
[C]/[0] = 1.08 in WASP-12b (our high [C]/[ 0] model case), the [O] abu ndance in the disk 
must have been depleted by ~ 50%. Recently, iMadhusudhan et al.l (l2011bl ) performed a more 
detailed analysis of the formation of WASP-12b and concluded that the primordial disk was 



deplete d in [O] by 41%. The discrepancy in our numbers arises because IMadhusudhan et al. 



(l2011bl ) used elemental abundances of the host star WASP- 12 (iFossati et al.ll2010bl ). which 
are 3.54 x 10"^ and 7.94 x 10"^ for [C] & [O], respectively. 

It is possible that the depletion of [O] in WASP-12b occured because the host star 
accreted fractionated refractory materials (that trapped 41% of [O], in the WASP-12 case) 



^It should be noted that, because the base level of the water clouds was not determined, the water 
abundance in the deep well-mixed regions of Jupiter is still unknown. However, this does not change the 
observed result that the mixing ratio of H2O gradually increases with depth. 
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from the protoplanetary disk during planetary formation. In our Solar System, planetary 
migration could have affected giant planet composition to some extent, as Jupiter and Sat- 
urn, in particular, are thought to have moved around considerably during planetary accre- 
tion JXsiganis et allboosl : IComes et aPboOsI : iMorbidelli et al.lboosi IWalsh et allboilj ). with 
Jupiter perhaps coming as close in as 1.5 AU. But, as far as we know, giant planets never 
migrated through the terrestrial planet region of our system. By contrast, in the WASP- 12 
system, planetary migration was evidently much more pronounced. WASP- 12b migrated 
from the outer parts of the nebula to its present location close to the star. Meanwhile, rocky 
planets formed in the hot, inner parts of the nebula may have migrated in even closer and 
been consumed by the star. If WASP-12b accreted additional material during its journey, 
that material would have be en depleted in O relative to C, possibly accounting for the high 
[C]/[0] ratio of the planet. iMadhusudhan et al. J2011bh mention this possibility, but they 
rule it out because they argue this would require that [C]/[H] in the envelope of WASP- 12b 
should be close to the host star's value, which is not the case. We don't agree with this 
objection, however, because all known planets accrete elements heavier than He more effi- 
ciently than they do H. Jupiter, fo r example, is enriched in C, N, and S compared to the 
Sun by a factor of 3 (lBeattylll999l ). Further observations may be needed to determine the 
validity of this mechanism. 

The accretion of refractory elements onto a star has been proposed as one of the rea- 
sons why solar twins and analogs in the solar neighborhood have enhanced h eavy elemental 
abundances compared to the Sun (IMelendez et al.ll2009l : iRamirez et al.ll2009l ). These studies 
found that the abundances of heavy elements in these solar analog stars increase with their 
condensation temperature. The authors attribute the apparent depletion of refractory ele- 
ments in the Sun to the existence of terrestrial planets, and they suggest that Sun accreted 
refrac tory-depleted mat erial from the nebula during the formation of the solar syste nfl. In- 
deed, [Chambers I (120101 ) showed that adding 4 of Earth-like and carbonaceous-chondrite- 
like material to the solar convection zone brings the Sun's elemental abundance in line with 
the mean abundances of solar twins. 



Recently, ISchuler et al.l ( 120111 ) derived precise elemental abundances for ten stars using 
high- resolution spectroscopy. All ten of their stars have at least one giant planet around them 
at different orbital distances. Their analysis indicated that four stars, all of which have hot 
Jupiters (~ 0.05 AU), show positive correlations between refractory elemental abundance 
and condensation temperature. This implies that these stars may have accreted refractory- 
rich planetary material or cores. If a similar accretion happened during the formation of 



^Note that the accreted material is not depleted in elements such as carbon because they have low 
condensation temperatures and so did not condense in the inner part of the solar nebula. 
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the WASP-12 system, then "pollution" signatures in the atmosp here of WASP- 12 may be 



obser vable. Formation models of protostars from molecular clouds (jWuchterl fc Tscharnuter 



20031 ) indicate that solar-mass protostars have thin convective envelopes (~ 0.02 Mq), similar 



to the present day Sun, and hence mixing of deposited material may not be significa nt. But 
it has been suggested that WASP-12b may be losing mass to its star (ILi et al.ll2010l ). If this 
is the case, then it could be difficult to determi ne how much of the r efractory material has 
accreted onto WASP-12 during its formation. iFossati et all (l2010bl ) performed a detailed 
spectropolarimetric analysis of WASP-12 to look for pollution signatures due to the material 
lost by WASP-12b. They found hints of pollution but were unable to draw firm conclusions. 
A differential analysis of WASP-12 twins (with the same effective temperature, age and 
metallicity), identifying their abundances with high precision, is required to determine if 
the refractory elemental abundance of WASP-12 does indeed increases with condensation 
temperature. 



5. Conclusions 

In this study, we analyzed how a disequilibrium mechanism such as photochemistry can 
affect the observed abundances of H2O, CO, CO2 and CII4 in the WASP-12b atmosphere. 
We considered two models, with [C]/[0] = 0.54 (solar) and [C]/[0] = 1.08 (twice solar). 
Although our photochemical results agree that high [C]/[0] is needed to explain the observed 
high abundance of CH4 and lack of observable H2O, they also indicate 2h2 and HCN are 
more abundant than CH4 and should be taken into consideration in modeling hot-jupiter 
atmosphere^. More importantly, our results indicate that 2h2 is the dominant absorber at 
1.6 and 8 micron in WASP-12b's atmosphere and the absorption features may possibly be 
arising from 2h2 rather than CH4. Observations with Hubble Space Telescope's WFC3 can 
resolve this discrepancy. 

We also propose a possible mechanism for the origin of the excess carbon observed in 
WASP-12b. Following other authors, we suggest that WASP-12 may have accreted rocky, O- 
rich material from the nebula during the formation of the system, leaving the disk relatively 
enriched in other volatiles such as C and N. WASP-12b then accreted some of this high 
[C]/[0] material, which thereby gave rise to the high [C]/[0] ratio of the planet. Testing this 
hypothesis requires that we understand whether WASP-12 is currently stealing mass from 
its planet, WASP-12b. A high precision abundance analysis of WASP-12 twins and analogs 



^"^Notc that this prediction is not spccificany a rcsuh of our photochemical model, as our thermodynamic 
equilibrium models predict this, as well. 
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can shed light on the refractory elemental abundance of WASP-12 and the possible origin of 
excess carbon in WASP-12b. 

The authors would like to thank the referee, Julianne Moses, for pointing to us the 
importance of 2h2 and HCN chemistry and for her in-depth analysis of our work which helped 
in improved photochemical models and our current manuscript. R. K and J.F.K gratefully 
acknowledge funding from NASA Astrobiology Institute's Virtual Planetary Laboratory lead 
team, supported by NASA under cooperative agreement NNH05ZDA001C, and the Penn 
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Fig. 1. — Temperature profile (red solid curve) and eddy diffusion profile (green dashed line) 
us ed in our photochern i cal m odel. The temperature profile is taken from one of the models 
of iMadhusudhan et al.l ( 120111 ) with no inversion, as a profile with inversion is ruled out by 
the data. This profile is then recalculated using hydrostatic equilibrium to be consistent 
with our photoch emical model. Eddy diffusion profile is taken from the dayside profile of 
Line et all J2010[ ). 
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Fig. 2. — Comparison of the normalized flux of a GOV stellar spectrum from iPickled (jl998[ ) 
stellar flux library (green curve), which is assumed for the star WASP-12, with an F2V star 
(top blue curve) and the Sun (red bottom curve). Inset figure shows the UV region. 




Volume Mixing ratio 

Fig. 3. — Equilibrium (dashed) and photochemical (solid) mixing ratio profiles of major 
species, for [C]/[0] = 0.54 (the solar value). Filled squares represent equilibrium values at 
the lower boundary. The mixing ratios refer to volume mixing ratio. The Helium abundance 
is 0.07836. 



1e-08 




Volume Mixing ratio 

Fig. 4. — Mixing ratio profiles from equilibrium (dashed) and photochemical (solid) models 
for [C]/[0] = 1.08 (twice solar) 




Fig. 5. — Approximate optical depths of the CH4 7.7 micron band (red sohd) and the 7.5 
micron 2h2 band (blue dashed) as a function of pressure. The optical depth of 2h2 is larger 
than CH4 indicating that it may be the major absorber in WASP-12b's atmosphere. 




Fig. 6. — Diagrams showing the rapid shift in species' concentrations as C/0 increases from 
0.99 (top panel) to 1.0 (bottom panel). 




Fig. 7.— Photolysis rates of H2O and CH4 for [C]/[0] = 0.54 (dashed, solar) and [C]/[0] = 
1.08 (sohd, super-solar) including 2/i2. 
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Fig. 8. — Photochemical mixing ratio profiles of H2O, CO, CH4 and CO? derived from our 
model of the daysid e atmosphere of HD 189733b. For comparison, see lLine et al.l (120 lOf ) and 
Moses et al.l toilh . 
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Table 1. Reaction list and rate constants fjZalinlell201ll ) used in this study. Only the 
forward rate constants are given as we calculate reverse rate constants from the forward 
rate assuming thermodynamic equilibrium (see discussion in Section 2). For three body 
reactions, the first and second row represent the low and high pressure rate limits, 
respectively. A full version of the table is available in the electronic edition of the 

Astrophysical Journal. 



Number 


React ants 




Products 


Rate*" 




Reference 


1 


H + H + M 


— ^ 


H2+ M 


8.8 X 10~^^(T/298)~ 


0.60 


Baulch et al.(1992) 




H + H 




H2 


1.0 X 10"^^ 






2 


+ H + M 




OH + M 


4.3 X 10-32 




Tsang & Hampson (1986) 




+ H 




OH 


1.0 X 10-12 






3 


H2 + 




OH + H 


3.5 X 10-i3(T/298)2-67e- 


-3160/T 


Baulch et al.(1992) 


4 


H+ OH + M 




H2O + M 


6.6 X 10-32(T/298)- 


-2.1 


Javoy et al.(2003) 




H+ OH 




H2O 


2.7 X lO-iOe-^5/^ 




Cobos & Troe (1985) 


5 


H2+ OH 




H2O + H 


1.6 X 10~^^{T/298y-^°e- 


-1660/T 


Baulch et al.(1992) 



2-body reaction rates are in cm^ s ^; 3-body rates are in cm^ s ^. 
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